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1. Introduction 

The charmonium spectrum consists of eight narrow states below the open charm threshold 
(3.73 GeV) and several tens of states above that. Below the threshold almost all states are well 
established, on the other hand very little is known at higher energy. In addition to that there have 
been discoveries [1] of several new charmonium-like states above the threshold. While some of 
them appear to be consistent with conventional charmonium others don't seem to behave like stan- 
dard meson and could be made of a larger number of constituent quarks. 
In the next sections we will review recent BaBar results on this field. 

2. Study of the 7/ y/(0 in two-photon interactions 

The Y(3940) was observed for the first time by Belle [2] in B decays and then confirmed by 
BaBar [3]. In a re-analysis [4] of the BaBar data sample the precision of the Y(3940) parameters 
was improved and was also found evidence for the decay X(3872) — > J/yco. This confirmed an 
earlier Belle claim [5] for the existence of this decay mode. A subsequent Belle paper [6] reports 
the evidence of a structure in the process yy—tJ/ yco that they dubbed the X(3915) with mass and 
width values similar to those obtained for the Y(3940) by BaBar [3]. BaBar has recently performed 
a study of the process yy — > J/yco [7] to search for the X(3915) and the X(3872) using a data 
sample of 5 19 fb _1 . Figure 1 shows the reconstructed J/yco mass distribution after all the selection 
criteria have been applied. A large peak at near 3915 MeV/c 2 is observed with a significance of 7.6 
a. The measured resonance's parameters obtained from a maximum likelihood fit are m x ^ 9l5 ^ = 
(3919.4 ± 2.2 ± \.6)MeV/c 2 , r x(3915) = (13 ± 6 ± 3)MeV. The measured value of the two-photon 
width times the branching fraction, is r rr (X(3915)) x ^(X(3915) ->■ J/ya) = 52 ± 10 ±3 eV 
and 10.5 ± 1.9 ±0.6 eV for the spin hypothesis J=0 and J=2, respectively, where the first error is 
statistical and the second is systematic. In addition a Bayesian upper limit (UL) at 90 % confidence 
level (CL) is obtained for the X(3872): r yr (X(3872)) x .^(X(3872) -> J/yco) < 1.7 eV, assuming 
J=2. 




m(JA|jci)) (GeV/c 2 ) 

Figure 1: The efficiency-corrected invariant mass distribution for the J/yco final state. The solid line 
represents the total fit function. The dashed line is the background conttibution. The solid histogram is the 
non J/yco background estimated from sidebands. The vertical dashed line is placed at the X(3872) nominal 
mass. 
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3. Search for the Z\ (4050)+ and Z 2 (4250)+ 

In 2008 the Belle Collaboration reported the observation of a resonance-like structure called 
the Z(4430) + decaying to ^f(2S)n + studying of the process B — > \y(2S)Kn [8]. This claim gen- 
erated a great deal of interest [9, 10] since such a state must have a minimum quark content 
ccdu, and thus would represent an unequivocal manifestation of four-quark meson state. The 
BaBar collaboration searched the Z(4430) + in an analysis similar to that of Belle in the pro- 
cess B — > y(2S)Kii and also in B — > J j\\fK% [11] but without finding any structure neither in 
the \j/(2S)n nor in the J /yn. Recently Belle performed an amplitude analysis of the J/ij/Kn [12] 
system finding no significant evidence of the Z(4430) + in agreement with BaBar. In 2009 the 
Belle Collaboration reported also the observation of two additional resonance-like structures simi- 
lar to the Z(4430) + in the study of the B° — > % c iK~7i + [13]. These new structures were labeled as 
the Zi(4050)+ and Z 2 (4250)+, both decaying to Xd^ + - 

BaBar , using a data sample of 429 fb -1 , has recent searched for the Z\ (4050) + and Z2(4250) + in 
the process 8° ->■ % C \K~K + and in the decay B + — > K^% c \iz + [14], with the % c \ — > J/yy. In the 
BaBar analysis the % c \ K + mass distribution, after background subtraction and efficiency-correction, 
has been modeled using the angular information from the Kn mass distribution that has been rep- 
resented using the Legendre polynomial moments. The excellent description of the % c \ft + mass 
distribution given by this analysis approach shows that there is no need for any additional res- 
onance to model the distribution. Figure 3 shows the result of the fit done on the Xci^ + mass 
spectrum using two or one scalar Breit-Wigners with parameters fixed to the Belle measurement. 
In all the fit cases there are no significant resonant structures, since the statistical significance ob- 
tained is very low, less than 2a. The Upper limits (ULs) at the 90 % CL on the branching fractions 
are for the one resonance fit: m{B° Z + K~) x m(Z + -»• Xc\K+) < 4.7 x 10~ 5 while for the two 
resonances fit: &(B° -»• Z+K~) x ^(Z+ -»• %c\X + ) < 1.8 x 10~ 5 and m{B° -> Z+K~) x m{Z+ -> 
Xc\X + ) <4.0x 10~ 5 . 

4. Study of the J/\{f7t + 7t~ system via Initial State Radiation (ISR) 

In 2005 BABAR discovered the Y(4260) in the process e + e~ ->■ yi SR Y (4260), with the 7(4260) - 
J/\j/7i + 7l~ [15]. Since it is produced directly in e + e~ annihilation it has J PC = 1 . The obser- 
vation of the decay mode J/xj/n^n [16] established that it has zero isospin. However it is not 
observed to decay to D*D* [17], nor to D*D* [18], so that its properties do not lend itself to a 
simple charmonium interpretation, and its nature is still unclear. A subsequent Belle analysis [19] 
of the same final state suggested also the existence of an additional resonance around 4. 1 GeV/c 2 
that they dubbed the Y(4008). BABAR has performed recently a new analysis [20] of this process 
using 454 fb . In this new study the region below 4.0 GeV/c 2 has been studied for the first 
time. As shown on Figure 4(a) in that region there is an excess of events above the J j\\f side- 
bands background. To understand the nature of this contribution a detailed study of the \\f(2S) line 
shape has been performed and the result shows that it is not possible to discount the presence of 
a //i//7T + 7r~ continuum cross section in this region. Figure 4(a) shows the fit to the J/\j/n + 7i^ 
mass distribution. A clear signal for the Y(4260) is seen; the values obtained from an unbinned- 
maximum-likelihood fit are: m Y ^ 2 60) = 4244 ±5 ±4 MeV/c 2 , rY (426 o) = 114+}^ ±7 MeV and 
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Figure 2: (a),(b) Background-subtracted and efficiency-corrected %c\K mass distribution for B — > % C \K%. 
(a) Fit with the Zi(4050)+ and Z 2 (4250)+ resonances, (b) Fit with only the Zi(4050)+ resonance. (c),(d) 
Efficiency-corrected and background-subtracted %ci ft mass distribution in the Kit mass region where Belle 
found the maximum resonance activity: 1.0 < m 2 (Kn) < 1.75GeV 2 /c 4 . (c) Fit with Z\ (4050) + and 
Z2(4250) + resonances, (d) Fit with only the Z(4150) + resonance. The dot-dashed curves indicate the fitted 
resonant contributions. 



T ee x 0§{Jly% + %-) = 9.2±0.8 ±0.7 eV. There is no evidence for the Y(4008) seen by Belle [19]. 
In the new BaBar analysis a detailed study of the 7T + 7r~ system from the Y(4260) decay to J / \j/n + 7l' 
has been performed. The n + 7i~ mass distribution shown in Figure 4(b) seems to peak around the 
/o(980) mass; however the peak is displaced from the nominal /o(980) position, since it is around 
940 MeV/c 2 . The fact that the peak is displaced and the particular shape of m(7t + 7t~) distribution 
may suggest a possible interference between the /o(980) and m(n + 7l~) continuum. To test this 
possibility the /o(980) line shape is taken from the BaBar analysis [21] of the — > 7i + 7t~7t + 
and this amplitude has been used in a simple model to describe the 7i + n~ mass distribution: 
\\fpol + £ !< ^/ (980)| 2 where pol is a polynomial function used to describe the m(7t + 7t~) contin- 
uum and ^/ (980) i s tne amplitude from the D$ — > [21] analysis; allows for a phase 
difference between these amplitudes. The result of this study is shown in Figure 4(b) and indicates 
that if there is a real /b(980) contribution to the decay of the Y(4260) to J j\\f% + %~ its contribution 
is not dominant: ^^Mm^z^JEl = (i 7 ± 13)%. 

5. Study of the \if(2S)7t + x~ system via Initial State Radiation (ISR) 

In addition to the Y(4260), two more J PC = 1 states, the Y(4360) and the Y(4660) have 
been reported in ISR production e + e~ -> ^(25 , )7r + 7T _ [22, 23]. While the Y(4360) was discov- 
ered by BaBar [22] and then confirmed by Belle [23] the Y(4660) was only observed by the Belle 
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m(J/i(f7C*ll-)(GeV/c 2 ) m(Jc + 7c)(GeV/c 2 ) 

Figure 3: (a): the J mass spectrum from 3.74 GeV/c 2 to 5.5 GeV/c 2 ; the points represent the data 

and the shaded histogram is the background from the J/y/ sideband; the solid curve represent the fit result, 
(b) the n + n~ distribution from the Y(4260) decay to J/\j/7Z + 7Z~ . The solid curve represent the fit using the 
model described in the text. 
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Figure 4: The y/(2S)7T + 7T~ invariant mass distribution from 3.95 to 5.95 GeV/c 2 for the y/(2S) -> 
J/y/K + n~; the points with error bars represent the data in the ty(2S) signal region, and the shaded his- 
togram is the background estimated from the yf(2S) sideband regions. The solid curve show the result of the 
fit. 



Collaboration. BaBar performed a new analysis [24] using all its available dataset collected at 
the T(nS), n=2,3,4; that corresponds to an integrated luminosity of 520fb~ 1 . The y(2S)n + K~ 
mass spectrum for the i//'(25) — > //i//7r + 7r~ is reported in Figure 4. BaBar observes two reso- 
nant structures, that have been interpreted as the Y(4360) and the Y(4660), respectively. The 
parameters values obtained from an unbinned-maximum-likelihood are for the first resonance 
m y(4360) =4340± 16±9 MeV/c , ry(436o) = 94 ±32 ± 13 Me V, and for the second onemy( 4660 ) = 
4669 ±21 ±3 MeV /c 2 , r y(4660) = 104 ± 48 ± 10 MeV. 
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